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Abstract 

Different families of generalized CS for one-dimensional systems with general time- 
dependent quadratic Hamiltonian are constructed. In principle, all known CS of sys¬ 
tems with quadratic Hamiltonian are members of these families. Some of the con¬ 
structed generalized CS are close enough to the well-known due to Schrodinger and 
Glauber CS of a harmonic oscillator, we call them simply CS. However, even among 
these CS there exist different families of complete sets of CS. These families differ 
by values of standard deviations at the initial time instant. According to the values of 
these initial standard deviations one can identify some of the families with semiclassical 
CS. We discuss properties of the constructed CS, in particular, completeness relations, 
minimization of uncertainty relations and so on. As a unknown application of the 
general construction, we consider different CS of an oscillator with a time dependent 
frequency. 
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1 Introduction 

1.1 General 

Coherent states (CS) play an important role in modern quantum theory as states that 
provide a natural relation between quantum mechanical and classical descriptions. They 
have a number of useful properties and as a consequence a wide range of applications, 
e.g. in semiclassical description of quantum systems, in quantization theory, in condensed 
matter physics, in radiation theory, in quantum computations, in loop quantum gravity, 
and so on, see, e.g. Refs. mmm- Despite the fact that there exist a great number of 
publications devoted to constructing CS of different systems, an universal definition of CS 
and a constructive scheme of their constructing for arbitrary physical system is not known. 
However, it seems that for systems with quadratic Hamiltonians there exist at present a 
common point of view on this problenQ. Starting the works [JJ; (U [D, ID] CS are defined as 
eigenvectors of some annihilation operators that are at the same time integrals of motion, 
see also [TTHTM3HT1IT1HT6]. Of course such defined CS have to satisfy the corresponding 
Schrodinger equation. In the frame of such a definition one can in principle construct CS 
for a general quadratic system. This construction is based on solutions of some classical 
equations, their analysis represent a nontrivial part of the CS construction. 

In this article we, following, the integral of motion method, construct different families of 
generalized CS for one-dimensional systems with general time-dependent quadratic Hamil¬ 
tonian. Analyzing these families, we see that some of them are more close to the well-known 
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due to Schrodinger and Glauber CS (see El) of a harmonic oscillator, we call them simply 
CS. However, among the latter CS there exist still different families of complete sets of CS. 
These families differ by values of standard deviations at the initial time instant. According 
to the values of these initial standard deviations one can identify some of the families with 
semiclassical CS, as was demonstrated by us in the free particle case m- We discuss proper¬ 
ties of the constructed CS, in particular, completeness relations, minimization of uncertainty 
relations and so on. As an application of the general construction, we consider CS of an 
oscillator with a time dependent frequency. 


1.2 Basic equations 

Consider quantum motion of a one-dimensional system with the generalized coordinate x 
on the whole real axis, x £ R = (— 00 , 00 ), supposing that the corresponding quantum 
Hamiltonian H x is given by a quadratic form of the operator x and the momentum operator 
Vx = - ihd x , 

H x = rip 2 x + r 2 x 2 + r 3 (xp x + p x x) + r 4 x + r 5 p x + r 6 , (1) 

where r s = r s (t), s = 1,..., 6 are some given functions of the time t. We suppose that these 
functions are real and both H x and p x are self-adjoint on their natural domains Dh x and 
D Px respectively, see e.g. [19] . 

Quantum states of the system under consideration are described by a wave function 
T (x, t) which satisfies the Schrodinger equation 

ihdt'b (x, t) = H x <5 ( x , t ). (2) 


In what follows, we restrict ourselves by a physically reasonable case r 4 ( t ) > 0. In this case, 
we introduce dimensionless variables, a coordinate q and a time r as follows 


,-1 I'* ds 2H l 2 

q = x, ~’ T = lm = rl n{s)ds ' T{t) = 2FMfy 


(3) 


where l is an arbitrary constant of the dimension of the length. The new momentum operator 
p and the new wave function if) ( q , r) read 


l ( ml 2 \ 

P=^Px = -id q , i/j(q,T) = Vm (lq,—Tj 


(4) 


so that I 1 ! (x,t)\ 2 dx = (< 7 ,t)| 2 dq. 

In the new variables, equation |2]) takes the form 

Si/j(q,T) = 0, S = id T — H, (5) 


where the new Hamiltonian reads 

p 2 

H = — + aq 2 + /? (qp + pq) + gq + vp + £ . 


( 6 ) 


Here a = a (r), /3 = f3 (r), g 




a ^ 2h 2 n (t) 
l r 5 ( t ) 


14 r 2 it) or _, s l 2 r 3 (t) 

^ T) = 2hWWY 


_ l 3 r 4 (f) 

2 fi2 n(t)’ 


^ - 'bw ^ 12 r e{t) 

" (T)= 2 hWWY £(r)= 2^MQ’ 


(7) 


are dimensionless real functions on r if t is expressed via r by the help of eqs. (j3j). In what 
follows, we call S the equation operator. 
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2 Constructing time-dependent generalized CS 


2.1 Integrals of motion linear in canonical operators q and p 

First we construct an integral of motion A (r) linear in q and p. The general form of such 
an integral of motion reads 

A (t) = / (r) q + ig (r) p + ip (r), ( 8 ) 

where / (r), g (r) and p (t) are some complex functions on r. The operator A(t) is an 
integral of motion if it commutes with equation operator (0) , 


S,A(t) 


= 0 . 


(9) 


In the case if the Hamiltonian is self-adjoint, the adjoint operator A 1 (r) is also an integral 
of motion, i.e., 


(r) 


= 0 . 


( 10 ) 


The commutator 


A(t) ,i f (t) 


reads 


A (r), (r) 


= S = 2 Re [g* (r) / (r)]. 


( 11 ) 


Substituting representation © into ©, we obtain the following equations for the func¬ 
tions f(r),g (t) , and ip (r): 

/ (r) + 2/3 (r) / (r) - 2ia (r) 5 (r) = 0, 
g (r) - if (r) - 2/3 (r) g (r) = 0, 

p (t) + 1/ (t) / (r) - ip (t) 5 (r) = 0. (12) 


It is enough to find the functions / (r) and g (t) , then the function (r) can be found by 
a simple integration. In addition, without loss of the generality we can set ip (0) = 0. 

Equations ED imply that S is a real integral of motion, S = const. In what follows we 
suppose that (5=1, which means 


Re \g* (r) / (r)] = Re [g* (0) / (0)] = 1/2. (13) 

Any nontrivial solution of two first equations ED consists of two nonzero functions / (r) 
and g (r). That is why we can chose arbitrary integration constants in these equations as 

/ (0) = c\ = |cr| e iM1 , g(0) = c 2 = |c 2 |e^, |c 2 | ^ 0 , |d| ^ 0 . (14) 

In terms of the introduced constants, condition ED yields 


M |ci|cos(/ii - H 2 ) = 1/2. (15) 

Under the choice (5=1, operators A(t) and d/ (r) become annihilation and creation 
operators, 

i(r),it(r)] = 1 . (16) 

It follows from eq. © and (TCH) that 

q = g*{ T ) A (t) — <p (r) + g (t) d f (r) — ip* (r) , 

ip = f* (r) \A (r) - p (r)j - / (r) [i f (r) - (r)l . (17) 
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I. We note that the two first equations can be reduced to a one second-order dif¬ 
ferential equation for the function g (r), such an equation has the form of the oscillator 
equation with a time-dependent frequency w 2 (r), 

9 (r) + w 2 (r) g (r) = 0, J 2 (r) = 2a - 4/3 2 - 2/3. (18) 

If we have an exact solution g (r) for a given function ui 2 (r), then the function / (r) can be 
found via the function g (r) as 


/ (t) = 2*/3 (r) g (t) - ig (r). 

One can chose the functions a (r) and /3 (r) such that 

w 2 (r) = 2a (r) — 4/3 2 (r) — 2/3 (r). 

For example, if we chose 

a ( T ) = 7/w 2 (t) > £ = 0 = v = e = 0, 

then we are dealing with Hamiltonian of the form 


H = 


P 


W 2 (r) ?2 


2 2 

II. In addition, the one-dimensional Schrodinger equation 

- d 2 4- (g) + y (g) 4- (g) = ^ (g), 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


can be identified with eq. (HU if 9 t, V (g) g (r), V {q) - E ^ w 2 (r). 

III. It should be also noted that two first equations m can be identified with a particular 
form of the so-called spin equation, see [20], 


iV = (<tF) V, V = 


(24) 


with 


F (r) = -- (2a + 1, i (2a - 1), 4i/3). 


2.2 Time-dependent generalized CS 

Let us consider eigenvectors | z,t) of the annihilation operator A(t) corresponding to the 
eigenvalue z, 

A(t)\z,t) = z\z,t) . (25) 

In the general case 2 is a complex number. 

It follows from eqs. (HU and (l25ll that 


q (r) = (z, t |g| z, r) = g* (r) \z - (r)] + g (t) [z* - ip* (r)] , 

ip (t) = (z, t \p\ z, t) = f* (r) [z-( p (r)] - / (r) [z* - p>* (r)] , 

- = / (t) 9 (r) + ig (t) p (t) + ip (r). ( 26 ) 


Using dm, one can easily verify that the functions q (r) and p (t) satisfy the Hamilton 
equations 


4(t) 


3H 
~d p ’ 


p{r) 


dH 

~dq’ 
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where H = H (q,p) is the classical Hamiltonian that corresponds to the quantum Hamilto¬ 
nian ©• Thus, the pair q (r) and p (r) represents a classical trajectory in the phase space 
of the system under consideration. All such trajectories can be parametrized by the initial 
data, qo = q (0) and po = p (0). 

Being written in the ^-representation, eq. reads 

[/ (t) q + g (t) d q + if (r)] (q\z, r) = z{q \z, t) . (27) 


General solution of this equation has the form 


(<?l z , t ) 


®z 1C2 {V, t) = exp 


f(r) g 2 

9(t) 2 


g( T ) 


q + x(r) , 


(28) 


where \ (r) is an arbitrary function on r. 

One can see that the functions >f>^ lC2 (q, t) can be written in terms of the mean values 
q (r) and p (r) given by eqs. 


K lC2 (?, t) = exp |ip (r) q - [q - q (r)] 2 + % (r) 

where x ( T ) is again an arbitrary function on r. 

The functions <f>° lC2 satisfy the following equation 


(q,r)=A(r) 4>r 2 (q,r) 


where 


A (r) = id T x (r) + aq 2 (r) - 


1 


p 2 (t) 


/ 


— ivp (t) — /3 ie. 


(29) 


(30) 


(31) 


If we wish the functions (1551) satisfies the Schrodinger equation (O, we have to fix \ ( r ) 
from the condition A (r) = 0. Thus, we obtain for the function % (r) the following result: 


X ( T ) = </H T ) + In N, 

<t> ( T ) = / I ( T ) ~ 


— ivp (r) — /3 — ie > dr 


p 2 (r) 


were N is a normalization constant, which we suppose to be real. 

The probability densities generated by the wave functions (1291) have the form 


(32) 


pT 2 (q,r) = |$r 2 (q,r) f = N 2 exp {- 


[q-q ( r )] 


2 Iff ( r )l 

Considering the normalization integral, we find the constant N, 

exp(—Re0(r)) 


2 Re (f> (r) 


pT 2 (q, t) dq = 1 => N = 


y/V^\g(T)\ 


(33) 


(34) 


Thus, normalized solutions of the Schrodinger equation that at the same time are eigen¬ 
functions of the annihilation operator A (r) have the form 


$T 2 (q,r) = 


1 


: exp < ip (r) q — 


\]\Z2tt \g (r)| 
and the corresponding probability densities read 

1 


PT 2 (9,r) = 


V^7t|5(t)| 


exp < 


/(r) [g-g(r )] 2 

9( T ) 2 


[q-q O )] 2 

2 \g ( r )| 2 


i Im <f> (7 


(35) 


(36) 


In what follows we call the solutions (1551) the time-dependent generalized CS. 
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3 Time-dependent CS of quadratic systems 

Using eqs. (ED and (1251) we can calculate standard deviations <r q (r), a p (r), and the quantity 
oqp (t), in the generalized CS, 

v q ( T ) = ({q- (g» 2 ) = \J (g 2 ) - (q) 2 = \g ( T )\, 

Vp 0 ) = \/((p - (p » 2 > = \j ( p 2 ) - ( p ) 2 = 1 / 0 ) 1 . 

°qp 0) = \ <(g ~ (g)) (p - (p» + (p- (p)) (g ^ (g»> 

= * [!/2 - g (t) /* (r)]. (37) 

One can easily see that the generalized CS (1351) minimize the Robertson-Schrodinger 
uncertainty relation ED, 

Vq 0) CTp - cr 2 qp (r) = 1/4. (38) 

This means that the generalized CS are squeezed states [3]. 

Let us analyze the Heisenberg uncertainty relation in the generalized CS taking into 
account restriction ED, 

0) cr p (r) | 2 Re ( c . C2 ) = ^ v^l + 4 (Im {gf*)) 2 > ^. (39) 

Then using (l37l) . we find a q (0) = a q = |c 2 1 and a v (0) = a v = |ci| , such that at r = 0 
this relation reads 


( 7 9 ( 7 pl2Re(cJ C2 ) = Y J + tMMsin^ 2 -^ 1 )] 2 . (40) 

Taking into account eqs. ED, we see that if jLti = = At the left hand side of (15(71) is 

minimal, such that 

GqGp — 1/2, Oqp = 0. (41) 

One can see that the constant fi does not enter CS ED. Then, in what follows we 
consider generalized CS with the restriction fii = fj , 2 = /i = 0. Namely, such states we call 
simply CS. 

Now restriction ED takes the form ci = |ci|, C 2 = | C 21 , 2 ci = c 2 1 , such that 

5 ( 0 ) = |c 2 | = a q , /(0) = |ci| =a p = 2—. (42) 

2(7 q 

Thus, a qp = a qp (0) = i [1/2 — 5 (0) / (0)] = 0, which is consistent with eqs. (HT1) . 

With account taken of eqs. ED, ED and ED, we obtain the following expression for 
the CS: 


$0 (g,T) = 

<ht)= r 

Jo 



exp < ip (r) q — 


ZO) [g-gO)f 

9 0) 2 


*Im^>(r) 


V 2 0 ) 


ZO) 


sO) 


— ivp (r) — /3 — ie > dr. 


(43) 


In fact, we have a family of CS parametrized by one real parameter-the initial standard 
deviation cr 9 > 0. Each set of CS in the family has its specific initial standard deviations 
u q . Different CS from a family with a given a q have different quantum numbers z, which 
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are in one to one correspondence with trajectory initial data go and pq. It follows from eq. 
that 

go „ Im z 


z = 


2a, 


+ ia q p 0: go = 2a q Re 2 , p 0 = 


The probability density that corresponds to the CS (l43l) reads 


Pz q {q, t) = 


1 


V2ira q (7 


■ exp - 


k-q(T)Y 

(r) 


(44) 


(45) 


One can prove that for any fixed a q states (l43l) form an over complete set of functions 
with the following orthogonality and completeness relations 


J C? (g, t)^” (g, r) dq = exp (^z'*z - 


Izf + RI 2 ' 


, Vr, 


(g, r) ^>z q (gb r)d 2 z = 7rS (q — q '), d 2 z = dRez dlmz, Vr. 


(46) 


4 An exact solution of oscillator equation with time- 
dependent frequency and related CS 

Let us consider the following function or (r), 

9/ 1 2 

u 2 (t)=uj 2 + , 2 ° , uj 2 <uj 2 (r)<uj 2 max , W (±oo)= W 2 , (47) 

cosh loqt 

where ui and loq are some positive constants, w max > uj. The function 0 J 2 (r) is an even 
function, which decreases monotonically as |r| changes from 0 to 00 , 

to 2 (± 00 ) = w 2 , uj 2 (0) = w^ ax =UJ 2 + 2ul . (48) 

For uj = 1 and ujq = 2 -1 / 2 , the plot of the function uj (r) has the form 

oj(t) 



The general solution of equation 1151) can be written as 
iA loq tanh(u}QT) 


SO) = 


B 




- Wo 


cos(wr) + 


UJ 1 + Wn 


— B uj 0 tanh(wor) 


sin(ccr) 


UJ 


(49) 
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The restriction (fl3l) , (fl4l) , and (1421) that set the CS from the entire set of generalized CS 
lead to the following relations for the constants A and B: 

B = g($) = \g (0)| = <j q , A = f(0) = |/ (0)| = l/2cr g . (50) 

Using eqs. fliil) and (l50l) we calculate the mean trajectories q(r) and p (t) according eqs. 

( ED ), 

q(T) = 2Re[g(r)z*] = g(T)\ A=A> B=B , 

p( T ) = g (t) | a=Ai b=b = q[r) , A = -cr ” 1 Im z, B = 2 a q Re 2 . (51) 

For w m in > 0 the mean trajectory q(r) can be presented as 

q(r) = RqR(t) sin[wr + 0(r) + 0 O ] , (52) 

where functions R(t) and 0(r) and constants Rq and 0 q are 


R(t) = \J 1 + tanh 2 uj 0 t , 1 < R(t) < w 1 y / a 


,2 . 
°0 i 


0(r) = arctan 


w 0 


tanh wqt 


Rq — 


\/pgw 2 + gg(w 2 + wg) 4 

(w 2 +wg ) 2 


—A < 0(r) < A, A = arctan ( 

V u) 


go 


sin©o = —, cos©o = 
Rq 


p 0 UJ 




Rq((jJ 2 + Wq) 


2\2 ' 


(53) 


Thus, we deal with a quasiharnronic motion with the frequency w and an amplitude that is 
changing in time in finite limits and with a time-dependent phase that is slowly changing in 
also finite limits. 

Let us derive the case of a harmonic oscillator with a fixed frequency io > 0 from the 
above formulas. To this end we have to set wo = 0 and a = w 2 /2, v = /3 = e = 0 such that 
w 2 (r) = w 2 . Then 


. . isinwr , . cos lot 

g (r) = <jq coswr + —-, / (r) = —-b t^usinwr, 

ZCTqLd ZCTq 

q(r) = qo cos lot + — sinwr, p(r) = po cos lot — uiqo sinwr, (54) 

u> 

and 2 = cr p qo + icr q po . Taking all that into account in Eq. (1431) . we obtain the following 
representation for CS (in the above given definition) of the harmonic oscillator: 


K q (g, ■ 



exp 


f 1 / (r) 

Z 

{ 2 3 (r) 

r / ( r ) J 


/* (r) z 2 
f(r) 2 



(55) 


For these CS 


(t) =a q ^j 1 + ^ 4fT 4 sin2 UT ’ a P ( T ) = a P\J l ~ (! “ 4a g w2 ) s in 2 wt, 


and the corresponding probability density reads 

1 


Pz q (g, t) = 


V^TTO-g (r) 


exp < -- 


[ g ~ g ( T )] 2 

2 <7 o (r) 


(56) 


(57) 


One has to consider the following three cases: 























a) <j q y/ 2w = 1, then 

<Tq (r) = <J q , CT p (t) = (Tp, <J q (t) (Jp (t) = O q (Jp = 1 / 2 , Vr. 

b) o q \f2u) < 1, then 

CT 9 ( r )lmin = CT 9 ( T )It=^ = ^, CT 9 Mlmax = ^9 

a P ( T )lmi„ = a P ( T )lr=2iti^ = CT 9 W > °P ( T )l m ax = G V ( T )I t =^ = °i» 

^W^WImin = CT 9 l T = 3 ffl = V 2 , 


CT 9 (t) (r)| max = ct 9 (r) <7, (r)l = 2n+i „ = 


c) a q \/2uj > 1, then 


1 + 4<j4o; 2 

8(JqU! 


-, n e N = 0,1,2,... . 


M T )Ln = = 


2<J q UJ 


CT 9 ( T )lmax = CT 9 WI T =Jff = ff 9 


°P ( T )lmin = a P ( T )lr=^ = CT P> CT P ( T )lmax = a P ( T ) l r= 2z!±I = ^9 W > 

(Tq (r) CTp (r)| min = CTg (r) ct p (t)| t= „ ^ = 1/2, 

1 + 4ct^w 2 

CT 9 (r)CT 9 (r)| ma x= CTg(r)CT, (r)| T= 2 j^i 5 = 8g2 ^ , nSN . 

We can see that in the case a) the Heisenberg uncertainty relation is minimized in 
the CS (l55l) . In the same case these CS coincide (up to a phase factor) with the well- 
known Schrodinger-Glauber CS nzi- CS with <jq obeying b) and c) minimize the Heisenberg 
uncertainty relation periodically, but the product a q (r) cr q (r) is always restricted by the 

limits 1/2 and 

1 8ct„oj 


Setting uio=<jj = a = h l = f3 = E = 0 1 and taking into account the limits 

iA 

lim g (t) = B cos (cut) -sin (ojt) , lim g (r) = B + iAr , 

UJq —^0 (jJ UJq jLO — ^0 


2 tanh(woT) 


lim c(t) = {iA- 

UJq 


B, 


we obtain from (|43l) CS of a free particle studied by us in the Ref. [18 / 
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